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ABSTRACT 

A  common  problem  found  in  advance  structural  materials  is  water  entrapment.  This  problem  is  a 
major  cause  of  material  degradation.  In  metals  it  can  lead  to  corrosion  and  in  composites  it  adds 
unnecessary  weight  to  the  structure  and  can  lead  to  material  degradation  especially  after  fi-eezing  and 
thawing.  Thermography  has  been  investigated  as  a  means  of  detecting  water  entrapment  in  metallic 
structures.  A  simple  model  has  been  derived  that  accurately  describes  the  thermal  response  of  these 
structures  to  short  heat  pulses.  In  this  paper  results  on  Aluminum  panels  with  various  amounts  of  water 
entrapment  will  be  presented.  Sensitivity  relations  will  be  derived  and  validated. 

Keywords:  NDE,  thermography,  thermal  models,  calorimetric  models,  defects,  water  entrapment 

L  INTRODUCTION 

The  detection  of  entrapped  water  behind  metallic  structures  is  a  difficult  problem  to  solve  in  a 
nondestructive  manner.  X-ray  radiography  is  commonly  used  to  detect  the  entrapped  water  but  the 
technique  is  not  very  portable  in  general  and  requires  special  safety  precautions.  Ultrasonic  attenuation 
methods  have  been  successfully  used  but  suffer  from  a  lack  of  lateral  resolution.  In  this  paper 
experimental  thermographic  results  and  simple  thermographic  models  will  be  presented  and  discussed. 

Simple  thermographic  models  have  been  previously  introduced  [1,2]  to  successfully  model  the 
surface  temperature  evolution  of  panels  with  flat  bottom  holes.  In  those  papers  the  hole  diameter  and  depth 
were  varied  to  simulate  different  amounts  of  mass  loss  due  to  corrosion.  In  all  cases  those  simple 
thermographic  models  sufficed  to  accurately  predict  the  maximum  contrast  temperatures  observed  in  the 
experiments.  This  paper  will  closely  follow  the  format  and  spirit  of  the  previous  references.  Once  again  it 
will  be  seen  that  simple  physical  arguments  go  a  long  way  in  describing  complicated  thermal  phenomena. 
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Fig.  I :  (top)  Drawing  indicating  key  parameters  used  in 
study.  Straws  of  various  diameters  (ranging  from  D  =  I/8” 
to  I/2”)  were  glued  to  aluminum  panels  of  different 
thickness  “to”.  The  straws  were  filled  with  various  amounts 
of  water  “d”  (Right)  Shows  a  single  thermographic  frame 
some  time  after  a  heat  pulse.  The  dark  regions  in  the  image 
represent  cold  spots  due  to  the  presence  of  water  in  the  back 
of  the  panel. 


In  this  paper  a  simple  theoretical  model  that  contains  all  the  relevant  parameters  required  to 
describe  the  thermal  transient  process  of  aluminum  panels  with  water  entrapment  will  be  presented.  Some 
of  the  parameters  used  in  the  model  include  the  amount  of  entrapped  water,  the  diameter  of  the  region  with 
entrapped  water,  lateral  heat  flow  effects,  anisotropic  properties,  finite  thickness  effects  of  the  metal  plate 
and  of  the  water  column.  In  spite  of  the  simplicity  of  the  model,  it  will  be  seen  that  the  model  fits  the 
experimental  data  remarkably  well. 


2.  EXPERIMENTAL  METHOD 


The  panels  used  in  this  study  were  imaged  using  a  commercial  pulsed  IRNDT  system 
(EchoTherm®,  Thermal  Wave  Imaging,  Inc.).  Flash  heating  was  provided  by  2  linear  xenon  flashtubes 
with  5  msec  flash  duration,  each  powered  by  a  6.4  kJoule  capacitor  bank.  The  system  was  equipped  with  a 
256  X  256  pixel  InSb  focal  plane  array  camera  (Radiance  HS,  Raytheon  Amber)  operating  in  the  2-5 
micron  spectral  range.  Continuous  12-bit  data  was  acquired  at  a  120  Hz  frame  rate  for  6  seconds  after 
flash  heating. 

The  material  used  in  this  study  was  aluminum  7075-T6.  To  model  water  entrapment,  straws  of 
different  diameter  were  glued  to  the  back  of  the  panel  and  filled  with  water  to  various  heights.  Fig.  1  (Left) 
shows  a  drawing  with  some  of  the  key  parameters  used  in  the  model.  The  parameters  “D”  and  “d” 
represent  the  thickness  of  the  panel,  the  diameter  of  the  straws  and  the  height  of  water  added  to  the  straws 
respectively.  Three  different  thicknesses  were  chosen  for  the  aluminum  panels.  Four  different  straw 
diameters  were  studied.  The  straws  were  filled  with  water  to  different  heights.  A  total  of  36  (3x3x4) 
different  conditions  were  studied.  Table  1  summarizes  all  the  parameters  used  in  fabricating  these  panels. 
The  front  side  of  the  panels  were  painted  with  a  flat  black  paint  to  avoid  emissivity  variation  problems. 


Fig.  2:  (Left)  This  fig.  shows  the  data  curves  require  to  generate  a  single  contrast  curve.  The  curve  labeled  “Water” 
represents  the  thermal  time  history  of  a  point  on  the  surface  of  the  aluminum  panel  directly  beneath  a  straw  with  filled 
with  water.  The  curve  labeled  “Dry”  represents  the  thermal  time  history  of  a  point  away  from  any  water.  The 


contrast  curve  wtis  generated  by  subtracting  the  previous  two  curves.  (Right)  This  figure  shows  three  contrast  curves 
for  three  identical  water  drops  but  for  three  different  aluminum  plate  thickness. 
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Fig.  3:  (Left)  This  fig.  shows  four  contrast  curves  for  a  given  aluminum  plate  thickness  and  straw  diameter  but  for  four 
different  water  heights.  (Right)  This  fig.  shows  the  relation  between  the  “maximum  temperature  contrast”  for  three 
aluminum  plate  thicknesses  and  four  straw  diameters. 

Fig.  1  (Right)  shows  an  actual  frame  taken  soon  after  the  thermal  flash  was  shot.  The  dark  areas  in  that 
photo  indicate  the  presence  of  water  in  the  back  side  of  the  panel.  The  straws  were  staggered  so  as  to  keep 
the  water  as  far  as  possible  from  one  another  in  order  to  minimize  interaction  effects  while  maximizing  the 
number  straws  in  one  panel. 


3.  EXPERIMENTAL  RESULTS 

In  a  standard  experiment,  digital  data  acquisition  begins  a  few  frames  before  the  capacitor  banks 
are  discharged  through  the  Xenon  arc  lamps.  In  our  experiments  the  frame  rate  was  set  to  100  frames/sec. 
A  total  of  660  frames  were  acquired  and  the  entire  experiment  lasted  over  6.6  sec.  Fig.  2  (Left)  shows  the 
entire  thermal  history  of  two  points  on  the  surface  of  the  panel.  The  curve  labeled  “water”  was  taken  from 
a  point  directly  below  and  in  the  center  of  a  straw  containing  some  amount  of  water.  This  curve 
characterizes  the  thermal  evolution  a  typical  site  that  has  some  water  over  it.  The  curve  labeled  “dry“  was 
taken  from  a  point  far  away  from  any  straw  containing  water  and  characterizes  the  thermal  history  of  a  dry 
aluminum  block.  Notice  that  the  thermal  history  curve  labeled  “water”  decays  faster  and  to  a  lower  value 
than  the  thermal  history  curve  labeled  “dry”.  This  is  due  to  the  fact  that  the  “water”  acts  like  a  thermal 
sink,  thereby  producing  a  decrease  in  the  surface  temperature.  Also  notice  that  the  “water”  curve  and  the 
“dry”  curve  have  a  region  at  early  times  with  very  elevated  effective  temperature.  This  region  of  the  plot  is 
difficult  to  model  because  it  includes  two  effects.  The  first  effect  is  the  direct  radiation  from  the  rapid 
temperature  rise  in  the  surface  of  the  test  panel  due  to  the  initial  flash.  The  second  effect  is  the  reflected 
radiation  from  the  walls  of  the  shroud  into  the  camera.  This  effects  disappears  when  calculating  the 
contrast  curve. 


The  difference  of  the  “water”  curve  from  the  “dry”  curve  is  termed  the  “thermal  contrast”  curve 
(shown  in  Fig.  2  Left).  Thermal  contrast  curves  start  and  end  with  zero  temperature  since  the  initial  surface 
temperature  and  final  equilibrium  temperatures  are  uniform  throughout  the  entire  panel.  Fig.  2  (Right) 
shows  three  contrast  curves  for  three  identical  water  drops  behind  three  different  plate  thicknesses.  It  is 
clear  from  that  figure  that  the  thicker  the  aluminum  plate  is  the  smaller  the  peak  thermal  contrast.  Fig.  3 
(Left)  shows  four  contrast  curves  for  four  water  drops  behind  a  1/32”  thick  aluminum  panel.  The  diameter 


water 
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Fig.  4:  This  Fig.  shows  a  schematic  representation  of  the  zero  lateral  flow  assumption  used  in  the 
model. 


of  the  straws  was  kept  constant  an  equal  to  ’/i”  while  the  amount  of  water  was  increased.  All  four  curves 
are  very  similar  in  shape  and  the  maximum  contrast  value.  When  the  experiment  was  repeated  but  with 
thicker  aluminum  plate,  the  four  contrast  curves  were  very  similar  again.  These  results  suggest  that  there  is 
some  form  of  saturation  effect  when  the  parameters  used  for  the  experiments  lay  in  the  range  as  shown  in 
Table  2.  It  is  important  to  note  that  an  ideal  contrast  curve  is  impossible  to  achieve  since  it  requires  two 
points  on  the  surface  of  the  panel  infinitely  apart.  Even  in  a  more  practical  sense  it  is  very  difficult  to 
generate  systematic  contrast  curves,  because  of  difficulties  in  selecting  equivalent  reference  point  on  the 
surface  of  the  panel.  Some  of  the  differences  seen  in  Fig.  3  (Left)  are  simply  the  result  of  difficulties  in 
finding  equivalent  points  over  all  areas  of  interest. 

Finally  Fig.  3  (Right)  shows  a  plot  of  the  maximum  thermal  contrast  as  a  function  of  the  straw 
diameter  for  three  aluminum  plate  thickness.  In  all  ceises  studied  the  effect  of  water  content  (or  water 
height  in  this  case)  was  negligible.  Each  data  point  on  this  graph  is  the  average  to  the  peak  contrast  for  the 
four  water  contents  studied  (d=l/l6”,  1/8”,  3/16”  and  1/4”).  The  error  bars  represent  two  standard 
deviations  of  all  peak  contrast.  In  some  case  the  standard  deviations  are  smaller  than  the  spot  used  to 
represent  the  data  point.  This  again  illustrates  the  fact  that  for  the  range  of  parameters  used  in  this  study, 
the  effect  of  water  height  was  insignificant.  We  will  see  later  on  how  our  simple  model  is  able  to  account 
for  this  effect.  The  larger  the  panel  thickness  the  smaller  the  contrast.  For  the  thickest  aluminum  panel 
(to=l/8”)  only  one  data  point  could  be  obtained  due  to  the  small  value  of  the  maximum  contrast.  For  the 
panel  with  thickness  of  1/16”  three  data  points  were  obtained. 

4.  CALORIMETRIC  MODEL  (ZEROTH  ORDER  APPROXIMATION) 

A  simple  theoretical  model  (zeroth  order  approximation)  has  been  derived  and  it’s  based  on 
simple  calorimetric  arguments.  Wet  and  dry  regions  are  defined  in  this  model  and  it  is  assumed  that  no 
energy  flows  between  them.  Fig.  4  shows  a  schematic  representation  of  the  model,  where  the  two  regions 
(dry  and  wet)  have  been  physically  separated  to  stress  that  there  is  no  energy  flow  between  them.  In  the 
next  section  a  more  refined  model  will  be  derived  that  takes  into  account  lateral  heat  transfer  effects. 

By  using  simple  calorimetric  arguments  we  can  write  that  q,  =  p,(A,tJ  c,  T,  and  q2=  Pi(A2to)c,T2 
+  P2(A2d)c2T2  The  parameters  q,,  A„  t^,,  p„  and  c,  represent  the  energy  deposited  over  the  area  A,  of  an 
Aluminum  panel  with  thickness  t^  density  p,  and  the  specific  heat  Cj.  The  parameters  q2,  Aj,  d,  pj,  and  C2 
represent  the  energy  deposited  over  the  area  A2  of  a  composite  panel  made  from  an  aluminum  panel  of 
thickness  t^  and  a  layer  of  water  of  thickness  d,  water  density  P2  and  the  water  specific  heat  C2.  It  was 
assumed  in  this  model  that  the  initial  temperature  of  the  coupon  was  uniform  and  equal  to  zero  degrees. 
Finally,  if  it  is  assumed  that  the  energy  deposited  on  the  surface  of  the  sample  per  unit  area  is  constant,  i.e., 
q,/Aj  =  q2/A2  =  Q,  then  the  temperature  difference  (or  thermal  contrast)  between  both  blocks  Tj  -  Tj  =  AT 
will  be 


4T.-2-±f - - ]  (,) 

P.'i  >.  + 

This  equation  correctly  describes  some  of  the  features  of  pulsed  thermography  when  applied  to  this 
problem. 

1 .  The  contrast  (AT)  increases  linearly  with  the  amount  of  energy  deposited  (Q). 

2.  The  higher  the  specific  heat-density  of  the  substrate  (PiCjT)  the  smaller  the  peak  contrast  (AT  i) 

3.  As  the  water  content  decreases  (d  0)  the  contrast  vanishes  (AT  ->  0). 

4.  As  the  amount  of  water  grows  indefinitely  (d  00)  the  contrast  saturates  (AT  Q/p,c,to). 

It  is  interesting  to  note  that  this  simple  formula  shows  a  strong  dependence  of  the  contrast  temperature 
AT  on  the  amount  of  water  “d”  inside  the  straw.  As  was  indicated  repeatedly  in  the  previous  section,  this 
was  not  an  experimentally  observed  fact.  It  is  anticipated  the  lateral  heat  conduction  effect  (not  included  in 
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this  model)  is  responsible  for  the  insensitivity  of  the  contrast  temperature  AT  on  the  variable  “d”.  In  the 
next  section  a  simple  model  that  takes  into  account  lateral  heat  flow  effects  will  be  presented. 

5.  CALORIMETRIC  MODEL  (FIRST  ORDER  APPROXIMATION) 

In  this  section  the  previous  model  will  be  modified  to  allow  for  lateral  heat  flow  effects.  A  “poor 
man’s  finite  element  approximation”  (i.e.  three  elements  only)  will  be  used.  All  the  elements  will  be 
thermally  interconnected.  Fig.  5  shows  a  schematic  representation  of  this  model.  In  this  model  it  will  be 
assumed  that  the  in-plane  thermal  conductivity  and  the  out-of-plane  thermal  conductivity  are  different. 
This  will  produce  the  most  general  results. 

In  this  model  it  will  also  be  assumed  that  all  the  energy  of  the  heat  pulse  is  absorbed  the  aluminum 
plate  of  thickness  As  a  result  of  this  heating  process,  the  temperature  “T^”  of  the  aluminum  plate  can 

be  derived  from  Q  =  pj-Ci  p-To  where  it  is  assumed  that  the  initial  temperature  of  the  panel  was  zero.  The 
quantity  Kl  represents  the  lateral  thermal  conductance  while  K^n-  represents  the  effective  out  of  plane 
contact  thermal  conductance.  The  lateral  thermal  conductance  can  be  expressed  in  terms  of  the  lateral 
thermal  conductivity  of  aluminum  k,  (which  in  this  case,  since  aluminum  is  isotropic,  is  just  the  thermal 
conductivity  of  aluminum) 


K,=k,.A^/R 

Xcff  =  h  •  Aj 


Similarly,  the  contact  conductance  can  be  expressed  in  terms  of  an  effective  contact  conductivity  “h” 
between  Aluminum  and  water.  A,  and  Aj  are  shown  in  the  fig.  5  and  represent  the  surface  area  of  the  dry 
region  and  the  wet  region  respectively.  A2  can  be  written  as  Aj  =  where  R  is  the  radius  of  the  straw 
while  A,  will  be  assumed  to  tend  to  infinity  00).  A  is  not  shown  explicitly  in  the  figure  but 

represents  the  lateral  cross  sectional  area  and  can  be  expressed  as  A^  =  27tR-to.  The  set  of  differential 
equations  that  define  this  problem  is 

HT  A 

p,.A,.t„c,.^  =  k,-^(T,-T,) 
dT  A 

Where  T„  Tj,  and  Tj’  are  the  temperatures  of  the  different  blocks  as  shown  in  Fig.  5.  This  set  of  coupled 
differential  equations  can  be  written  in  matrix  form  as 


Fig.  5:  This  figure  shows  the  building  blocks  of  our  simple  model. 
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^  =  kRT  (5) 
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This  set  of  three  coupled  linear  differential  equations  can  be  solved  by  first  diagonalizing  the  matrbc  R  . 
The  second  step  is  to  transform  the  problem  into  the  eigen-space  where  the  equations  get  de-coupled.  In 
this  space  the  solutions  are  simple  decaying  exponential  functions,  the  exponent  being  related  to  the  eigen¬ 
values  of  the  matrix  R  .  Finally  the  solution  in  the  eigen-space  has  to  be  transformed  back  into  the  original 
coordinate  set.  After  performing  the  previous  steps  and  imposing  the  following  boundary  conditions 


T(t  =  0)=  T. 
.  0 


where  T„  is  the  initial  temperature  of  the  aluminum  block  (which  can  be  derived  from  Q  =  p,-c,  p  T„)  the 
solutions  can  be  obtained  and  the  contrast  yields 


PiC|  to  ■^|(l  +  a  +  Ty  -4ar 


k  A,  I  pet 

where  a  =  — — - —  and  r  =  ^ — -  (don’t  confuse  the  variable  “t  =  time”  with  the  parameter 

h  A2  R  P2C2  d 

\  =  panel  thickness”).  The  parameters  and  X2  are  the  non-trivial  eigen-values  of  the  matrix  R  defined 
above. 

The  maximum  or  peak  thermal  contrast  can  be  calculated  by  differentiating  Eq.  7  and  the  result  gives 


AT^=- 


1  +  a  +  r  -  J(1  +  a  +  r)^  -  4ar  2V(u.+r)*-4ar 


Pi^i  t^(l  +  a  +  r  +  y(l  +  a  +  r)  -4ar)  [l  +  a  +  r  +  -^(l  +  a  +  r)  -4ar 
which  happens  at  a  time  give  by 


t  -Ml 
tpek-  h 


t^  j  l  +  a  +  r  +  ^(l  +  a  +  r)^ -4ar 

(1  +  a  +  r)^  -4ar  I  l  +  a  +  r”y(l  +  a +  r)^  -*4ar 


It  can  be  shown  that  in  the  limit  when  the  lateral  heat  flow  effects  disappear,  i.e  kL  0  (or  a  0)  then  Eq. 
8  reduces  to  Eq.  1  as  it  would  be  expected.  In  practice  the  contact  conductivity  “h”  is  much  smaller  that 
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the  lateral  thermal  conductivity  of  aluminum  “k”,  as  a  result  a  »  r  >  1.  If  eqs.  7,  8  and  9  are  expanded 
around  “a”  in  this  limit,  it  can  be  shown  that 


AT(t)  = 


h  I 

- 1 

P2C2‘1 


-e 


PR  r2 


(10) 


and  finally 


^Tpcak 


_Q _ L 

PiC,  a  t. 


(11) 

(12) 


It  is  important  to  stress  that  the  previous  three  equations  were  derived  in  the  limit  a  »  r,  which  happens  to 
be  the  “material  parameter  space”  region  considered  in  this  work  as  shown  in  Table  L  As  an  example,  a 
vanishing  amount  of  water  in  a  straw,  i.e.  when  d  «  t„  would  not  fall  within  these  “material  parameter 
space”  region  and  therefore  would  not  be  properly  modeled  by  the  previous  equations.  A  desirable 
property  of  eq.  1 1  is  that  it  depends  very  weakly  on  the  variable  “d”.  This  weak  dependence  on  “d”  was 
experimentally  observed  and  was  described  in  section  3.  The  variable  “d”  only  enters  into  eq.  1 1  via  the 
parameter  “r”.  If  “r”  «  “a”  then  (r/a)'^‘  1  and  AT  actually  becomes  an  independent  function  of  “d”. 


6.  ANALYSIS  OF  RESULTS 


Figure  6  (left)  shows  all  experimental  data 
points  and  the  fit  to  them  using  eq.  11.  It  can  be  seen 
from  the  graph  that  this  simple  model  fits  the 
experimental  results  fairly  well.  The  largest  fitting 
discrepancies  happen  for  the  panel  with  thickness 
1/32”.  Figure  6  (right)  shows  the  straw  diameter 
data  curves  for  three  different  panel  thickness  and  a  fit 
to  them  using  eq.  13.  It  can  be  seen  again  that  the 
fitting  results  for  the  1/32”  thick  panel  show  the  largest  discrepancies.  The  starting  parameters  that  were 
used  for  all  these  fits  are  shown  in  table  1  and  were  obtain  from  standard  material  properties  handbooks 
except  for  the  contact  conductivity.  The  final  fitting  values  were  found  to  be  reasonably  close  to  these 
starting  parameters.  The  starting  value  for  the  contact  conductivity  was  chosen  to  be  the  thermal 
conductivity  of  water  per  unit  centimeter.  The  final  fitting  value  was  found  to  be  significantly  larger  than 
this  initial  estimate. 


Aluminum  7075  T6 

Water 

p,  =  2.8  g/cm^ 

P2  =  1  g/cm^ 

Cj  =  0.84  J/(g  K) 

c,  =  4.1  J/(gK) 

k,  =  1.29  W/(cm  K) 

h  =  8.2  mW/(cm^  K) 

Table  2:  Relevant  material  property  values 


Maximum  Contrast  and  Fit 


Fig.  6:  This  figure  shows  a  fit  (solid  lines)  using  eqs.  8  or  1 1  to  the  experimental 
contrast  data  for  three  different  plate  thicknesses.  The  error  bars  represent  the  range  in 
peak  contrast  temperatures  for  the  various  amounts  of  water  used  in  the  experiments. 


It  is  important  to  emphasize  some  of  the  problems  found  when  dealing  with  the  type  of  data 
described  in  this  paper  and  when  using  this  simple  model.  The  first  problem  that  one  is  faced  with  is  in 
generating  the  contrast  curves.  As  was  mentioned  in  section  three,  and  ideal  contrast  curve  should  be 
generated  by  subtracting  a  reference  curve  taken  far  away  from  any  region  of  the  panel  that  has  water  from 
a  curve  taken  under  a  region  with  water.  This  requirement  was  difficult  to  meet  because  of  the  finite  size 
of  our  sample  eind  because  of  the  many  regions  with  water  that  we  had  in  our  panel.  To  overcome  this 
difficulty  one  could  use  one  aluminum  panel  for  each  water  drop  geometry.  The  problem  with  this 
approach  would  be  that  a  large  number  of  panels  would  be  require  and  non-systematic  errors  would  be 
introduced  to  the  data  such  as  room  temperature  variations  and  varying  amounts  of  heat  deposited  onto  the 
surface  of  the  sample.  Another  source  of  random  errors  was  the  straw  fabrication.  Straws  were  not 
perfectly  circular  and  the  diameter  as  a  result  had  random  errors. 

The  model  used  in  this  work  makes  some  drastic  approximations  such  as;  1)  a  very  small  number  of 
blocks  with  average  temperatures  to  describe  the  transient  behavior  of  the  system,  2)  temperature 
independent  material  properties,  3)  no  cooling  effects  introduced  through  radiation,  4)  uniform  initial 
temperature  distribution,  5)  single  drop  on  panel  (as  opposed  the  multitude  of  drops  actually  used  in  the 
panels),  6)  transverse  conductivity  approximated  by  a  contact  conductance  and  7)  does  not  take  into 
account  the  thermal  properties  of  the  straw.  It  is  anticipated  that  many  refinements  to  the  model  can  be 
achieved  by  improving  on  some  of  the  approximations.  Despite  its  crudeness  the  model  is  robust  enough 
to  fit  the  experimental  data  accurately  enough. 

7.  CONCLUSION 

Two  simple  models  have  been  developed  (figs.  4  and  5)  that  to  a  first  order  approximation 
describes  the  main  features  of  thermal  pulse  analysis  when  applied  to  a  planar  flaw.  The  first  model,  zeroth 
order  approximation  does  not  take  into  account  lateral  heat  conduction  effects.  The  second  model  takes 
into  account  lateral  heat  conduction  effects,  thickness  effects,  water  size  effects,  density  effects,  thermal 
transients.  The  model  correctly  predicts  the  relationships  between  the  previous  parameters.  Eqs.  7,  8  and  9 
are  the  main  output  of  the  mode  or  their  expansions  around  he  parameter  “a”  (eqs.  10,  1 1  and  12)  .  This 
relations  were  shown  to  model  correctly  the  time  dependence  of  the  thermal  contrast,  the  peak  thermal 
contrast  and  the  time  at  which  the  thermal  contrast  peaks. 
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